Universal formula for the muon-induced neutron yield 
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The experimental data on the yield of muon-induced neutrons for liquid scintillator, iron and 
lead accumulated during 60 years of muon interaction underground study have been analyzed. A 
universal formula connecting the yield with muon energy loss in the matter and neutron production 
in hadronic and electromagnetic showers is presented. 
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Introduction. -In the last decade a renewed interest to 
the problem of yield Y n of muon-induced neutrons be- 
comes evident. This is due to both the increased re- 
quirements to the accuracy of background definition in 
underground experiments and the growth of computing 
resources. The yield dependence on both the mean muon 
energy E^ and atomic weight A of the medium has been 
investigating using the FLUKA and GEANT simulation 
packages and their versions [l| , [2| , [3| , [4( . The yield cal- 
culation results are significantly different (Fig.l in [2[). 
The normalization of calculations to the experimental 
data becomes more complicated due to large straggling 
in these data. Until now no expression exists for Y n 
that binds together the muon energy deposition, nuclear 
properties of the matter and the neutron production pro- 
cesses in hadronic (h) and electromagnetic (em) showers 
generated by muons and developing in the matter. To 
calculate a yield value the approximate empirical laws 
Y n = p E xE* (for fixed A) or Y n = paxA? (for fixed 
Efj,) are used. The constants a, f3 are defined based on 
the results of calculations. Numerical fitting coefficients 
Pe, Pa are entered to get the agreement between calcu- 
lations and a set of experimental data available. 

The form of dependence Y n {E tl ) = aE was proposed 
in [91 . As follows from results of the measurements [fj, 



0,11,0 and calculations 0, 0, 0, [ijj, [ll|, the 

a value is in the range from 0.7 to 0.9. The values of 
exponents a and j3 representing the contributions of the 
neutron production channels play an important role in 
the analysis of the experimental and calculated data. 
Experimentally, the yield Y n is given by 



Y n = =^-{n/g/cm 2 ), 
I lip 



(1) 



where N n is the number of neutrons produced by a 
muon at energy E^ on the path length l^ in the matter 
with density p. E^ and I M are mean values for muon 
flux at a depth H and muon track lengths, respectively. 
The yield is connected with the medium properties and 
the characteristics of reactions of neutron production by 
expression: 



where N Q is the Avogadro number, {vu) is a mean 
value of the product of photonuclear /zA-interaction 
cross section and neutron multiplicity v, A is the atomic 
weight of the medium. Equation $Z§ follows from the 
dependence N n on (z/<r) and 7 M 
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(3) 



where ca [cm 3 ] is a concentration of nuclei A. 

Experimental data. -The table lists the measured yield 
Y LS for liquid scintillator (LS), Fe (Y Fe ), Pb (Y Pb ). The 
table covers almost all the data obtained for about 60 - 
year measurements of the neutron yield in underground 
experiments. The data are listed in the order of increas- 
ing of energy E^ to which authors attributed their re- 
sult. An error in determining average muon energy E^ 
was only shown in Ref. [12j. To estimate the E^ value 
error in other experiments summarized in the table we 

have used the expression 5E^ = 2/ y-Ey. It covers both 

the uncertainty E^ calculations at different sets of pa- 
rameters offered in Ref. [13j, [14|,[15| and deviations 
of E^ values from the E^ (H) dependence which can be 
seen in the table. The Mote Carlo calculations carried 
out recently [16|, [17J have resulted in a revision some 
Y L 



values. The most part of the LS data [6j, [8|, [llj . 

[19|, [20j, 21( was obtained using scintillator C n H2n 
n = 9.6 p = 0.78g/cm 3 0, [ll|, [U, fUJ (in the table 
the refined value Y^s = 4.1xl0~ 4 n/p/(g/cm 2 ) from [22| 
is included). 

It should be noted that the paper [3|| cites wrongly 
(Table IV in 0) the results of Ref. 0, @], [13, namely, 
out of 15 values taken from these works and included 
in the Table IV, seven values do not correspond to the 
published original data. The correct values of H, E^ and 
Yls are presented in the table of given paper and also in 
Ref. 01. 
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TABLE I: Measured neutron yield 









Y n xW- 4 ,n/fi/(g/cm 2 ) 








E^GeV 


H, m.w.e. 


Yls 


Yf, 


Y Pb 


Ref. 


year 


10.0±6.3 a 


20 


- 


0.98±0.01 


2.43±0.13 


[26] 


1954 


10.0±6.3 a 


60 


- 


- 


4.8±0.6 


[27] 


1970 


11.0±6.6° 


40 


- 


1.32±0.30 


4.03±0.36 


[9] 


1971 


13.0±7.2 


20 


0.20±0.07 


- 


- 


[18] 


1995 


16.5±8.1 


32 


0.36±0.03 


- 


- 


[19] 


2000 


16.7±8.2 


25 


0.47±0.05 
0.36±0.05 6 


- 


- 


[6] 


1973 


17.8±8.4° 


80 


- 


1.69±0.30 


5.66±0.36 


[9] 


1971 


20±9" 


110 


- 


- 


6.8±0.9 


[27] 


1970 


40±12.6 a 


150 


- 


3.31±0.96 


11.56±1.1 


[28] 


1968 


86±18 


316 


1.21±0.12 
0.93±0.12 6 


- 


- 


[6] 


1973 


110±21 a 


800 


- 


- 


17.5±3.0 


[29] 


1970 


125±22 


570 


2.04±0.24 
1.57±0.24 6 


- 


- 


[20] 


1986 


260±8 


2700 


2.8±0.3 


- 


- 


[12] 


2010 


280±33 


4300 


- 


- 


116±44 


[30] 


1973 


280±33 


3100 


4.1±0.5 


16.4±2.3 


- 


[21] 


2005 


280±33 


3100 


3.2±0.2 


19.0±1.0 


- 


[11] 


2011 


385±39 


5200 


c n+0.95 
J -°-1.02 

4.1±0.6 6 


20.3±2.6 




[8] 


1989 



"Vertical flux 
'Corrected values 



r 



The measurements were carried out in a global muon 
flux at different depths and energies E^ from 16.7 GeV to 
385 GeV. The experiments [6] , [20] detected the neutrons 
produced only in the counter LS; the results of Ref. [8j, 



|ll| . [2l| covered the neutrons generated in LS and iron 



of the setup structures (LS and iron masses were almost 
equal). 

The counters were located close to the mine ceiling of 
gypsum (Sp = 16.7 GeV) or salt (E^= 86 GeV) in the 
experiment [6J and close to the mine ceiling of salt (E^ 
= 125 GeV) in Ref. H. As a result of the Monte Carlo 
calculations in Ref. 16], it was obtained that the contri- 
bution of neutrons produced by shower particles in the 
standard rock around the detecting volume LS (C12-H26) 
enlarges the measured yield Yls by ^30%. Taking into 
account this fact and disregarding the small difference 
between compositions of LS and rock in experiments [6| , 
[201 ] and calculations [16[ , we have obtained the corrected 
values of Yls which are presented in the table. 

The similar correction is not suited for the results of 



21], 22| since in these experiments inner counters 



are detecting the neutrons produced in an inner volume 
of setup consisting of LS and iron in the same proportion 
as the peripheral part of the setup. 

With LSD and LVD which are almost identical in 



design and detection technique the yield was measured 
under different conditions of neutron detection: a) with 
inner counters crossed by a muon (LSD [8]), b) with all 
counters of the inner setup volume crossed by a muon 
(LVD [21|, |22j), c) with inner counters fired by any 
trigger pulse, including the muon trigger (LVD 23]). 
Here a muon is meant both a single muon and a muon 
group with shower accompaniment or without it. In 
all papers using the LSD and LVD data the yield was 
defined by formula 



Y, 



N' 



del 



LS 



where N det = JV||* 



NfipLshsr] 



(4) 



Np* is the number of detected 



neutrons, including produced in LS (Nls) and iron 
(Np e ), while N det = N L sVls + N Fe r] Fe , where r] LS , r] Fe 
are corresponding neutron detection efficiencies; Q is the 
fraction of neutrons produced in LS, N^ is an amount 
of muons, Ils is a mean length of muon tracks in LS. 
N ,Nfi and Ils have been determined directly in the 
experiment. N^ has been determined with due regard 
to the multiplicity of muon groups established using the 
tracking system data 2l[ . The Q fraction was calculated 



with the assumption that n = 7]ls = rjp e . The values of 
Q, 0.61, 0.60, 0.85 and n = 0.60, 0.90, 0.60 have been 
used for cases a), b) c), respectively. Case c) leads to 
the selection of neutrons at energies above 10 MeV and, 
as a consequence, to significant reduction in Yls what 



was indicated in [2l|, [2J|. For this reason, the result of 
23] is not included in the table. 



The recent Monte 
showed that rjLS^VFe 
the formula (j4|): 



Carlo calculations in Ref. [17| 
This leads to the need to change 



Y, 



N 



del 



LS 



Q 



NppLShs Q 7 1LS + (1 - Q)VFe ' 

Given fraction Q the yield Yp e can be also defined 



(5) 



Y P 



N 



del 



1-0 



Nfj.pFelFe QVLS + (1 - Q)VFe ' 



(6) 



where iFe is a mean length of muon tracks in iron and 
PFe is the iron density. 

New Q values were calculated in [25J for cases a), b) of 
the neutron detection with LSD and LVD. The Q frac- 
tion depends on the ratios of masses ku — Mls/Mfc, 
surface areas ks = SFe/Sps calculated per counter, 
atomic weights Ua = Als/Afc, and the exponent j3 



k^kMks 



(1 



k^kMks) 



(7) 



On the basis of the LSD experimental data in [25] it 
was found that j3 = 0.95 and Y Fe = 20.3 xl0~ 4 . To 
do this, the data of inner counters of first LSD level 
were used. In addition, these counters were detecting 
the neutrons produced by muons in a 8-cm thick steel 
platform beneath the setup. Using the data of the inner 
LSD counters of the second level, the fraction Q = 0.138 
has been determined for case a). The Y^g value of 
(4.1±0.6)xl0~ 4 corresponds to this fraction at efficien- 
cies n LS = 0.45, T)p e = 0.10 and /3 = 0.95. In the LVD 
experiment the value Q = 0.18 and the corresponding 
yields Y LS = (3.3±0.5)xl0- 4 , Y Fe = (16.4±2.3)xl0- 4 
were obtained under detection conditions bL Efficiencies 
7]ls = 0.75, r/Fe = 0.65 were taken from [ljj]. Thus, the 
Yls values from the reviewed papers exceed the corrected 
magnitudes by ~ 30 % (the table, Fig. [T]). 



The recent LVD results have been presented in [111 ]: 
Y LS = (3.2±0.2)xlO- 4 ,V fe = (19±l)xl0~ 4 . The yield 
values were obtained based on the data of counters 
without triggering pulses to avoid some methodical ef- 
fect reducing neutron detection efficiency in a counter 




FIG. 1: Dependence of the neutron yield on muon energy for 



scintillator. The curve is a function Y n 



4.03xl(r 6 £ 



fitting the experimental points (filled circles), open stars are 
uncorected data. 



with triggering pulse. In this case the efficiencies are 
rj LS = 0.0075, rjpe = 0.0107. All quantities (Q, 77, 7) 
except for the starting number N det of the detected neu- 
trons were calculated by the Monte Carlo method. 

Formula for the muon-induced neutron yield.- The 
data collected in the table including early measurements 



m, m 



£9j, [30( with iron and lead were an- 
alyzed using the conventional approach: a and j3 are 
constants independent of E^ and A, correspondingly. 
Using independence of a on E^ , for any A we can 
reduce the yield Y(E^) values to a certain arbitrarily 
chosen energy E and calculate the average value of 

(y(E/0): for £* = 16.7 GeV (Y LS ) = 0.34 (averaged 
over nine values), (Yp e ) — 1.70 (averaged over seven 
values), (Ypb) — 6.33xl0~ 4 n/p/(g/cm 2 ) (averaged 
over eight values). The ratio (Yls) / '{Yp e } is consistent 
with fi = 0.95, while (Y LS )/(Ypb) with (3 = 0.97 and 
(Y Fe )/(Yp b ) with /3 = 1.00. The large /3 values in the 
last two cases are mostly associated with excessive yield 
Y Pb = 116xl0~ 4 in experiment 0. 
The table data presented in Fig. [2] can be described by 
the expression 



Y n (A,E tl ) = cA f % 



(8) 



where (3 — 0.95, c is constant. Using the independence 
of /3 on A and assuming (3 — 0.95, the Yp e (E fi ) and 
Ypb(E^) data sets can be reduced to the Yls(E ii ) set 
(Fig. [21 lower panel). Fitting the yield set of 24 values 
Y L s(En) by expression Y LS = c(10.3) 095 ^ we get the 
best agreement with the data at c = (4.4±0.3)xl0 -7 



and a — 0.78±0.02. The same values c and a, but 
at larger uncertainties, result from corrected LS data 
(Fig. [TJ nine values). 

The constant c is close to the value of the rel- 
ative muon energy loss in nuclear interactions 
bh = 4.0xl0~ 7 (g/cm 2 ) -1 . Therefore, c is a rela- 
tive muon energy loss for neutron production c = b n 
and it has the dimension (g/cm 2 )^ 1 . Since neutrons are 
produced mainly in the em- and h- showers the constant 
b n should be associated with electromagnetic energy 
loss of muAgafonova2.texons 6 7 (em-shower generation 
mostly by means of bremsstrahlung) and nuclear loss 
bh (generation of h-showers where the main part of 
neutrons is produced). In the range from ~ 100 GeV to 
the extreme average muon energy underground ~ 430 
GeV the generation of em- and h-showers is proportional 
to E a . The bh value does not depend on E u and 
weakly depends on A: bh = 4.0 xlO -7 for standard rock, 
4.2 xlO -7 for water [l3(; the 6 7 loss varies slightly from 
12.1xl0~ 7 up to 14.2xl0 -7 for rock and from 8.2xl0~ 7 
up to 9.9xl0 -7 for water [13(. The 6 7 loss depend on the 
medium as Z 2 /A. The proximity of the values b n and 
bh reflects, on one hand, the dominant role of h-showers 
in neutron production and, on the other, the practical 
constancy of b n in a wide range of E^ and A. The 
values of exponents a,/3 in equation © are determined 
by neutron production processes in showers: in em- 
showers Y n (xE^ [7|, in h-showers Y n ccE u [8[, [lOj . 
[20] . Therefore, the resultant values a = 0.78, /3 = 0.95 
and b n — 4.4xl0 -7 obtained above are associated with 
the contributions of all neutron production processes, 
namely the shower generation by muons and the neutron 
production in showers via irA, NA, 7 J 4-reactions. 
Since the product b n E" [GeV/ p/ (g / cm 2 )] defines the 
muon energy loss for the neutron production, then, due 
to the yield dimension [n/p/(g/cm 2 )], factor A@ has the 
dimension of [n/GeV]. 

The yield value is contained in the formula for the neu- 
tron production rate r n = I fl (H)pAY n (E a , A) (n/cm 3 c), 
where I^H) (p/cm 2 c) is a muon intensity at a depth 
H and pa is a medium density. Using this formula one 
can write the expression for rate R n of muon-induced 
neutrons in the detector and its shield consisting of 
different materials. The neutron rate for a material Ai 
of a volume Vi and mass m, is given by 



Rni = Vir n = I a (H)p Ai ViY ni = I u (H)miY ni (n/c). (9) 
For all materials of the detector and shield we have 




150 200 



300 350 'ICO 

Muon energy, GeV 



FIG. 2: Dependence of the neutron yield on A and E^. Up- 
per panel: experimental points for lead (A=207, open circles), 
iron (A=56, filled circles) and scintillator (A=10.3, open tri- 
angles); the curves are functions Y = cA^E^ at different 
A and c = 4.4xl0" 7 ,/3 = 0.95, a = 0.78. Lower panel: 
neutron yield for scintillator; the experimental data for iron 
and lead are reduced to scintillator, the curve is a function 
F LS = 4.4xl0- 7 10.3°- 95 £?' 78 



for the heavy material (Fe, Pb) can be used for experi- 
mental determination of E^ at any overburden topogra- 
phy and rock composition. The accuracy of the proce- 
dure mi ght be not worse than in finding E^ by formulae 
in Ref. [13j, 14|, [15j. An approximation with constant 
parameters b n — 4.4xl0~ 7 cra 2 /g, a — 0.78,13 — 0.95 al- 
lows to use the formula ([5]) to calculate the yield for any 



En and A in underground experiments. Since all nu- 
clear effects produced by muons in the matter, including 
the production of radionuclides, are proportional to the 
neutron yield value the formula (jHJ is universal. How- 
ever, the magnitudes of the parameters are determined 
by the contributions of nuclear and electromagnetic pro- 
cesses and therefore albeit weakly but depend on E a and 
A. Due to the increasing requirements to the accuracy 
of the background determination in underground exper- 
iments the study of the neutron yield remains actual. 
This work was supported in part by the Russian Foun- 
dation for Basic Research grants 12-02-00213-a, 12-02- 
12127-ofi-m and SSh-871. 2012.2. 
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R n = I a {H)^ mi Y nz = I n (H)b n E«Zm t AP(n/c). (10) 

As it follows from equation (jHJ the neutron yield is highly 
dependent on E u {ccE n ) and A((xA - 95 ). So, its value [3] 
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